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I investigate an SU(3) gauge model with 12 fundamental fermions. The physically interesting 
region of this strongly coupled system can be influenced by an ultraviolet fixed point due to lattice 
artifacts. 1 suggest to use a gauge action with an additional negative adjoint plaquette term that 
lessens this problem. 1 also introduce a new analysis method for the 2-lattice matching Monte Carlo 
renormalization group technique that significantly reduces finite volume effects. The combination 
of these two improvements allows me to measure the bare step scaling function in a region of the 
gauge coupling where it is clearly negative, indicating a positive renormalization group /3 function 
and infrared conformality. 



Gauge models with many fermions or fermions in 
higher representations can develop a conformal phase 
characterized by the emergence of an infrared fixed point 
(IRFP) in the gauge coupling. Both the conformal sys- 
tems and those that are still chirally broken but are 
very near to the conformal window could be relevant for 
physics beyond the standard model. During the last few 
years many of these models were studied using various 
lattice simulation techniques [TVH2]. The theory with 
SU(3) gauge fields and 12 flavors of fundamental fermions 
has been the subject of extensive investigations, but its 
infrared behavior is still controversial. Refs. [TJ |3] used 
the Schrodinger functional approach to numerically cal- 
culate the renormalization group /3 function and con- 
cluded that the theory has an IRFP, i.e. it is confor- 
mal. Ref.0 considered the system at finite temperature 
and reached similar conclusions. Both of these works 
used unimproved or only moderately improved actions 
at strong gauge couplings where lattice artifacts could 
seriously effect the results. At the same time studies of 
spectral quantities appeared to be more consistent with 
a chirally broken system[31[S]. Early studies using Monte 
Carlo renormalization group (MCRG) techniques were 
not able to push deep enough into the strong coupling 
and remained inconclusive [Sj. Recently a large scale 
study [El concluded that high precision data of spectral 
quantities prefer the chirally broken interpretation, but 
other groups interpret the same data as more compatible 
with the conformal behavior|13]. 

In this work I revisit the 12 flavor SU(3) system using 
MCRG methods. Due to two improvements, one in the 
lattice action, the other in analyzing the MCRG data, I 
am able to cover a wider coupling range and can demon- 
strate that in the investigated region the renormalization 
group (3 function (actually its lattice analogue, the bare 
step scaling function) has the opposite sign of an asymp- 
totically free theory, signaling the existence of an infrared 
fixed point and the conformality of the system. 

The basic observation that led to the modified action 
is the existence of an ultraviolet fixed point due to strong 
coupling lattice artifacts. It is well known that the pure 



gauge SU(iVc) theory both with — 2 and 3 exhibits 
a first order phase transition in the fundamental-adjoint 
plaquette gauge action space [HJ [TS] . This line ends in a 
second order point that has a (most likely trivial) ultra- 
violet fixed point (UVFP). For notational convenience 
I call this new fixed point UVFP-2, while I use G-FP 
to refer to the perturbative Gaussian fixed point at zero 
gauge coupling. While the first order phase transition 
and the UVFP-2 are lattice artifacts and independent 
of the G-FP and the continuum limit defined there, their 
existence can strongly influence, even completely change, 
the scaling behavior of the lattice model. 

Ref. \W[ studied the scaling of several observables of 
the SU(3) fundamental-adjoint pure gauge system with 
adjoint coupling Pa = 0, -2.0 and -4.0, far away from 
the endpoint of the first order line that occurs around 
/3a ~ 2.0. Nevertheless the data showed very large scal- 
ing violations, even lack of scaling, at couplings near the 
extension of the first order phase transition line. As ex- 
pected, the scaling violations decrease with negative ad- 
joint terms in the action, i.e. farther from the second 
order endpoint. Refs. [iTl fT8] studied the RG fiow lines 
in the pure gauge SU(2) system. They found that near 
the extension of the first order phase transition line along 
the fundamental plaquette action the RG flows away from 
the UVFP-2 and turns around sharply at negative adjoint 
coupling. This again indicates that in this region the sys- 
tem is strongly influenced by the FP associated with the 
second order phase transition. In a recent work [19 we 
studied the pure gauge fundamental-adjoint SU(2) sys- 
tem with the 2-lattice matching MCRG method. Our re- 
sults show that near the first order line and its extension 
towards negative adjoint couplings the MCRG matching 
method breaks down, the system is no longer in the basin 
of attraction of the perturbative G-FP. 

When two UVFPs exist, numerical simulations have to 
stay in the vicinity of either one of them to describe the 
corresponding continuum physics. Fortunately the basin 
of attraction of the UVFP-2 is at fairly strong coupling, 
and present day QCD lattice simulations are sufficiently 
far from it. This, however, might not be the case in many 
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FIG. 1: The approximate location of the phase transition 
/crossover in the fundamental-adjoint plane. The solid red 
line indicates first order phase transition while the dashed 
red line corresponds to crossovers. 

fermion systems where interesting physics is expected to 
occur at strong gauge coupling. 

Large number of fermions could change the phase 
structure of the pure gauge system so I started by the 
study of the phase diagram of the fundamental-adjoint 
plaquette gauge action with 12 fermions. I used nHYP 
smeared staggered fermions [20] and measured the pla- 
quette, the specific heat through the derivative of the 
plaquette, the Polyakov line and the chiral susceptibility 
on 8^ X 4, 8"*, 12^ x 4 and 12^ x 6 lattices. The spe- 
cific heat gave a very clear signal for a first order transi- 
tion, continuing along a crossover line, as indicated by the 
solid and dashed red lines in Figure [T] This phase transi- 
tion/crossover has no dependence on the temporal lattice 
size, it is a bulk feature of the system. In the crossover 
region the chiral susceptibility gives no signal at all and 
I found no evidence for a finite temperature phase tran- 
sition even at very strong gauge coupling. Otherwise the 
phase diagram of the 12 flavor system looks very similar 
to the pure gauge one. There is a first order line ending at 
a second order fixed point around {/3p 2 A, 13a ~ 3.6). 
At smaller (Sa there is a crossover that gets weaker with 
decreasing adjoint coupling. By Pa = —1.4 (the last 
point along the dashed red line) there is only a very weak 
signal left. I should note that I did not determine the 
phase transition with high precision - my goal was to 
establish the qualitative features of the phase diagram. 

The horizontal blue line a.t Pa — Figure [T] shows the 
region I studied in Ref. [S]. MCRG matching became 
impossible at stronger couplings, to the left of the blue 
line. In retrospect that was most likely due to the nearby 
crossover region. It is useful to recall that the leading 
order perturbative relation between the gauge coupling 
and the lattice couplings is 



This suggests that the coupling {pp, Pa) = {5.0,0.0) corre- 
sponds, at least perturbatively, to (/?_f, /3^)=(10.0, -2.5). 
The latter point is quite far from the crossover along the 
(3a/ Pf — —0.25 action line, indicated by the second blue 
line in Figure [T] Finally, the green line in the figure cor- 
responds to Pa/ Pf = —0.50, the limit where the adjoint 
plaquette overtakes the fundamental one and flips the 
system into a new universality class. 

If the basin of attraction of the perturbative G-FP is 
limited by the first order/crossover line, and Eq. [l]is any 
indication of constant physics, than along the Pa/ Pf — 
—0.25 line one could reach considerably stronger cou- 
plings than with the Pa — fundamental action. I have 
chosen this action for the investigation described in this 
paper. This is a rather arbitrary choice, and other ratios 
could work equally well or even better. 

The 2-lattice matching MCRG method is a powerful 
tool to numerically calculate the bare step scaling func- 
tion, the discretized lattice analogue of the RG P func- 
tion. The method has been used for many years and re- 
cently it has been discussed in detail in Refs.jnUHj. Here 
I summarize it only briefly, concentrating on the new de- 
velopments. 

The bare step scahng function is Sb{P;s) — P — P' , 
where P and /?' are gauge couplings with lattice correla- 
tion lengths related as ^(/3) = s^{P'), where s > 1 is an 
arbitrary scale parameter. In my work I always consider 
s — 2 transformations and in the following I will drop 
the index s in Sb- The lattice correlation length can be 
deflned through different observables, which leads to sys- 
tematical uncertainties in the definition of si,. These can 
be controlled by approaching the UVFP, just like in the 
case of the renormalized step scaling function. In this 
work I will not perform this continuum extrapolation, 
rather include the systematical uncertainties in the error 
estimate of the final result. The step scaling function ap- 
proaches a constant at the G-FP, vanishes at other fixed 
points (both UV and IR), and has the opposite sign of 
the RG P function where it is non-zero. 

The 2-lattice matching MCRG method relies on 
matching observables after several RG blocking steps. Its 
advantage is that the simulations do not have to be per- 
formed on volumes with lattice size comparable or larger 
than the correlation length, and most of the finite size ef- 
fects can be cancelled by comparing blocked observables 
measured on identical blocked volumes. 

Under repeated renormalization group blocking the ac- 
tion first flows towards the renormalized trajectory (RT), 
than along it, away from the UVFP. If the RG flow lines 
originating at P and /?' hit the RT at the same points 
but require Ub and rif, — 1 blocking steps to do so, the 
requirement on the correlation length, ^(P) — 2^(/3'), is 
satisfled. There are two steps to achieve this: 
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• Matching: If the blocked actions are identical, ob- 
servables measured on Ub times blocked conflgu- 
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FIG. 2: The optimized A/? at /3f = 6.5 with the action 
Pa/ Pf = —0.25 . The red circles are the finite volume cor- 
rected predictions, the blue diamonds have no finite volume 
correction in the optimization step. The left side of the fig- 
ure shows results after comparing blocking steps ni, = 3/2/1 
on different volumes, while right side is the result after one 
more blocking step. The data points are labeled by the final 
blocked lattice size L^. 

rations generated at /3 have to match observables 
measured on n;, — 1 times blocked configurations 
generated at (3' . Since each blocking step reduces 
the lattice size by a factor of two, the final measure- 
ments are usually done on very small volumes. In 
order to minimize finite size effects it is best to do 
the simulations on twice as large lattices at /3 than 
at /?' so the final measurements are done on the 
same blocked volume. Thus the shift in the gauge 
coupling is defined as A/3ci(/3; rtf,, Lh) = /3 — /?' if 

(0(/3;nb,ib)) = (0(/3';nfc-l,Lfc)), (2) 

where (O) denotes the expectation value of some 
short distance operator and Ljy is the volume after 
Ub and rifc — 1 blocking steps (the same for both 
sides of Eq. [2]). 

• Optimization: The quantity A/3 defined in the pre- 
vious step can differ significantly from the step scal- 
ing function s;, if the RG flow does not reach the 
RT in n;, — 1 steps. Most RG block transforma- 
tions have a free parameter, usually denoted by a, 
that can be optimized to minimize the number of 
RG steps needed to reach the RT. The optimized 
parameter is defined as the one where consecutive 
blocking steps predict the same shift, 

A/3ci(/3;»^6,i6,aopt) = A^o(/3;nb - l,Lf„aopt) • (3) 

To minimize finite size effects A/Jq on the two sides 
of Eq. [3] should be calculated on the same blocked 
lattice size Lb- Previous studies did not take this 
volume dependence into account and usually satis- 
fied Eq. [3] on different volumes. The error intro- 
duced this way is much smaller than the one intro- 



duced by not matching the volume in Eq. [2] but 
still it can be important when A/^o itself is small. 

Eq. |3] requires the comparison of simulations on 
three different volumes. It is easiest to illustrate this 
with a specific example. Let's assume we simulate on 
32^* volumes at some (3 value. After blocking the lat- 
tices nb = 3 times we measure observables on = 4 
lattices. We match these to observables measured on 
rif, = 2 times blocked lattices at some coupling and 
find A/3o(/3;n6 = 3, if, = 4) = /3 - P' . Since Lb = 4, 
simulations at /?' must have been done on IG'* lattices. 
Optimization requires that 

A/3o(/3; m = 3, it = 4) = A^o(/3; Ub = 2, Lb = 4) . 

To calculate the quantity on the right hand side we have 
to do simulations on 16'' and 8* volumes at /3 and (3' . 
Identifying the optimal RG transformation and corre- 
sponding A/3ci with Tifc = 3/2/1 blocking steps requires 
simulations on volumes 32^, 16'* and S''. The procedure 
can be repeated with different O operators and the stan- 
dard deviation between the predicted A/J^ values charac- 
terize the systematical errors of the matching. Results on 
larger volumes with more blocking levels provide further 
consistency checks. 

The rest of this paper illustrates the optimiza- 
tion/finite size correction process and shows the step scal- 
ing function for a range of gauge couplings. The simula- 
tions were done on volumes between 32* and 4* using 12 
nHYP smeared staggered fermions. The gauge action is a 
combination of fundamental and adjoint plaquette terms 
with fixed Pa/ Pf = —0.25 ratio. The lattice fermion 
masses were am = 0.0025 on the 32*, 0.005 on the 16* 
and 24*, 0.01 on the 12* and 8* and 0.02 on 6* and 4* 
volumes. The masses were chosen such that their val- 
ues match if they scale with their engineering dimension. 
This is not the right scaling if the anomalous mass hap- 
pens to be large. However these bare fermion mass values 
are so small that the data show no mass dependence even 
with masses twice as large as used here. For all practical 
purposes these mass values can be considered to be in the 
chiral limit. I used an RG block transformation based on 
2 HYP smearing steps with fixed inner parameters and 
considered 5 different operators as described in Ref.[i8| 

Figure [2] illustrates the optimization at Pf = 6.5. 
The left side of the figure shows the optimal A/3 after 
Ub = 3/2/1 blocking steps and final blocked lattices of 
Lb —2, 3 and 4. As is easy to check, the original vol- 
umes in the matching sequences were 32 — >■ 16 — >■ 8, 
24 ^ 12 6 and 16 -J' 8 -J> 4. The red circles show the 
results of the optimized matching, a consistent value be- 
tween all three volume series. The blue diamonds show 
the predicted A/3 without finite volume correction in the 
optimization. The result on the smallest volume set is 
clearly off, signaling large finite volume effects. The two 
larger volumes show very little deviation, it appears that 
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FIG. 3: The bare step scaling function. The different sym- 
bols correspond to predictions from optimized matching on 
different lattice volumes and blocking levels. 



and that made simulations possible at physically stronger 
gauge couplings. Second I corrected for a previously ig- 
nored finite volume effect that made the results obtained 
on different volumes after different blocking levels con- 
sistent. This finite volume correction also reduced the 
systematical errors that come from matching 5 different 
operators. The same approach could easily be applied to 
other models near the conformal window. 

The numerical calculations of this work were carried 
out on the HEP-TH computer cluster at the University 
of Colorado and at FNAL, under the auspices of USQCD 
and SciDAC. This research was partially supported by 
the US Department of Energy. 



at least with my blocking transformations and 5 opera- 
tors a finite volume of Lf, = 3 is already sufficient to min- 
imize these second order finite volume effects. The right 
side of the figure shows A/3 after one more blocking step, 
with Tih = 4/3/2. The largest volume in this case was 32* 
with the final blocked volume Li, — 2. Again, the finite 
volume corrected optimized data is significantly different 
from the uncorrected one but both are consistent with the 
Lb = 2 results of the left hand side. The error bars on the 
data points come from a combination of statistical and 
systematical errors. They are dominated by systematical 
errors in the nf, — 3/2/1 sequence and by statistical er- 
rors in the rtf, = 4/3/2 one. Comparing the finite volume 
corrected optimized results for A/3 on all three volume 
sequences and after n\j = 3/2/1 and = 4/3/2 blocking 
levels one finds Sf,(/3 = 6.5) = -0.15(2) . 

Figure|3]show A/3 (or Sb) at a range of gauge couplings. 
The red diamond points are from 16 — > 8 — > 4, ri;, = 
3/2/1, the blue crosses are from 32 16 ^ 8, rib ~ 
3/2/1, and the black circles are from 32 — > 16 ^ 8, 
Ub — 4/3/2 optimized matching. Where all three data 
points are available, they are consistent. Overall, the 
data show that Sb is negative in the investigated region, 
indicating that the RG /3 function has crossed zero and 
the measurements are on the strong coupling side of the 
IRFP. One should note that the data in Figure [3] does 
not correspond to any given RG transformation. Each 
point has a slightly different optimization parameter and 
can have a different IRFP as well. 

In summary. Figure [3] gives strong evidence that the 12 
fermion SU(3) system is infrared conformal. This is not 
the first MCRG investigation of this theory, but previous 
ones were inconclusive. The success this time had two 
sources. I considered an action farther away from a sec- 
ondary UVFP caused by strong coupling lattice artifacts 
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